Advanced workflow

Dynamic simulation control

ResSimAD.jl allows flexible control over simulation runs.

API functions ResSimAD.time_step and ResSimAD.step_to enable simulating the model for a single timestep, or to a specific time, respectively.

API functions ResSimAD.change_well_mode, ResSimAD.change_well_target, ResSimAD.shut_well and ResSimAD.add_well enable dynamic modification of well setups during a simulation run.

After changing well settings, time step size will be cut down to sim.scheduler.dt0. We can also use ResSimAD.change_dt to change the time step size.

Here is an example where well control settings are changed dynamically during simulation runs. This is very convenient in the context of well control optimizations.

# Create simulation model
using ResSimAD: get_model, change_well_target, change_dt, step_to, time_step

sim, options = get_model("example1");

# Change BHP control every 100 days
tsteps = [100., 200., 300.];
for t in tsteps
    p1_bhp = rand()*400. + 5500.;
    i1_bhp = rand()*400. + 6100.;
    change_well_target(sim, "P1", p1_bhp)
    change_well_target(sim, "I1", i1_bhp)
    change_dt(sim, 5.0)
    step_to(sim, t)
end

# Change BHP every 5 time steps
for i = 1:2
    p1_bhp = rand()*400. + 5500.;
    i1_bhp = rand()*400. + 6100.;
    change_well_target(sim, "P1", p1_bhp)
    change_well_target(sim, "I1", i1_bhp)
    change_dt(sim, 5.0)
    for j = 1:5
        time_step(sim)
    end
end

## Plot results
using ResSimAD: get_well_rates, get_state_map, get_grid_index, isproducer
using Plots
using Plots.PlotMeasures
cmap = cgrad(:jet);

# BHP
t = get_well_rates(sim, "P1", "TIME");

pw_p1 = get_well_rates(sim, "P1", "WBHP");
pw_i1 = get_well_rates(sim, "I1", "WBHP");
p1 = plot(t, pw_p1, color=:black, marker=true,
     xlabel="Day", ylabel="BHP (psi)",
     title="P1 BHP");
p2 = plot(t, pw_i1, color=:black, marker=true,
     xlabel="Day", ylabel="BHP (psi)",
     title="I1 BHP");
plot(p1, p2, layout=(1,2), legend=false, size=(720,280), bottom_margin = 10px)

# Oil rate and water injection rate
qo = get_well_rates(sim, "P1", "ORAT");
qw = get_well_rates(sim, "I1", "WRAT");
p1 = plot(t, qo, color=:black, marker=true,
     xlabel="Day", ylabel="Oil rate (STB/Day)",
     title="P1 oil rate");
p2 = plot(t, -qw, color=:black, marker=true,
     xlabel="Day", ylabel="Water inj. rate (STB/Day)",
     title="I1 water injection rate");
plot(p1, p2, layout=(1,2), legend=false, size=(720,280), bottom_margin = 10px)

# Define a function for plotting wells
function plot_wells_2d(plt, sim; markersize=5, color=:white)
    for w in values(sim.facility)
        i, j, _ = get_grid_index(sim.reservoir.grid, w.ind[1])
        if isproducer(w)
            marker = :circle
        else
            marker = :dtriangle
        end
        scatter!(plt, [j,], [i,], m=(marker, markersize, color), legend=false)
    end
end

# Plot state maps
po = get_state_map(sim, "po", t[end]);
sw = get_state_map(sim, "sw", t[end]);
p1 = heatmap(reshape(po, sim.nx, sim.ny), color=cmap, title="Po at day $(t[end])");
plot_wells_2d(p1, sim);
p2 = heatmap(reshape(sw, sim.nx, sim.ny), color=cmap, title="Sw at day $(t[end])");
plot_wells_2d(p2, sim);
plot(p1, p2, layout=(1,2), size=(700,450))

We can also add new well dynamically. This can be useful for field development optimization.

using ResSimAD: add_well

new_well = Dict([
    ("name", "I2"),
    ("perforation", [(8,12,1)]),
    ("radius", 0.5),
    ("mode", "bhp"),
    ("target", 6500.)
])

add_well(sim, "injector", new_well);

t_end = 1200.;
step_to(sim, t_end);

# Plot production curves
qo = get_well_rates(sim, "P1", "ORAT");
t2 = get_well_rates(sim, "I2", "TIME");
qw1 = get_well_rates(sim, "I1", "WRAT");
qw2 = get_well_rates(sim, "I2", "WRAT");
p1 = plot(t, qo, color=:black, marker=true, label = "P1",
     xlabel="Day", ylabel="Oil rate (STB/Day)",
     title="P1 oil rate");
p2 = plot(t, -qw1, color=:blue, marker=true, label="I1");
plot!(p2, t2, -qw2, color=:red, marker=true, label="I2",
     xlabel="Day", ylabel="Water inj. rate (STB/Day)",
     title="Water injection rate");
plot(p1, p2, layout=(1,2), size=(720,280), bottom_margin = 10px)